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Abstract. A novel technique based on Schwinger's proper time method is applied to the Casimir problem of 
the M.I.T. bag model. Calculations of the regularized vacuum energies of massless scalar and Dirac spinor 
fields confined to a static and spherical cavity are presented in a consistent manner. While our results agree 
partly with previous calculations based on asymptotic methods, the main advantage of our technique is 
that the numerical errors are under control. Interpreting the bag constant as a vacuum expectation value, 
we investigate potential cancellations of boundary divergences between the canonical energy and its bag 
constant counterpart in the fermionic case. It is found that such cancellations do not occur. 



PACS. 12.20.DS Specific calculations - 12.39. Ba Bag model 

1 Introduction tractive force was derived by Casimir as early as 1948 |Q. 

About two decades later, great efforts were made to inves- 
The effects of a classical static background field on the ob- tigate more complicated geometrical arrangements of the 
servables of a relativistic quantum field theory have been boundaries for a variety of free field theories and bound- 
investigated in detail for quite some time. Perhaps the a ry conditions ||. In particular, following the develop- 
easiest way of distorting a free quantum field is to impose ment of the bag models of hadrons |,|,|, there has been 
boundary conditions on a static surface. For the special increased interest in vacuum energies arising as a con- 
case of an electromagnetic field, obeying boundary condi- sequence of the boundary conditions on a static sphere, 
tions on two parallel, uncharged and static plates, an at- More recently, such calculations have been extended to 
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include interacting, renormalizable quantum field theories but in practice constraints on the computational effort can 

in the framework of perturbation theory 0. cause the errors to remain quite large. 

There are a number of approaches to calculating vac- i n the framework of the M.I.T. bag model there 

uum expectation values of field operators in relativistic is f or eac h ne i d ( sca lar, fcrmionic and vector) a linear 

quantum field theories. In general, one can classify these boundary condition that determines the eigenmodes and 

into two categories: local or global methods @. The global their energies. There is also a nonlinear boundary condi- 

mode summation techniques avoid the explicit occurrence t ion that guarantees the balance of the pressures at the 

of divergences by analytical continuation (Zeta function boundary. The latter is usually taken into account in an 

regularization §, heat kernel expansion §), whereas lo- average fashion, by minimizing the total energy of the 

cal methods take advantage of the fact that residual local bag with respect to the bag radius. Comparing the re- 

divergencesf] manifest themselves only as nonintegrable in- smt f this minimization with the experimental data, a 

finities of density functions on the boundary @. phenomcnological value for the bag constant is obtained. 

The main purpose of this paper is to present a consis- However, the nonlinear boundary conditions of the model 

tent calculation of the regularized vacuum energy of mass- can also be used to interpret the bag constant as a vacuum 

less scalar and Dirac spinor fields, that are confined to a expectation value which makes this additional and rather 

static and spherical cavity, using the local stress tensor artificial minimization condition redundant . 

method fBl XO . llfl - The coefficients of the leading bound- Tk • • • i ( u T TT 

ol — Pl_^ 6 I his paper is organized as follows: In section II, we 



ary divergences of the cavity energy have already been iUustrate the local stress tensor method in the case of the 

found elsewhere @,|l||3) based on asymptotic expan- straight . forward examp i c f sca i ar and Dirac spinor nc i ds 

sions. However, as one of these authors points out @, satisfying lincar boundary conditions on a plane. Well- 

the validity of this method can be questioned. It is there- knQwn resultg are reproduced using the momentum-space 

fore useful to check these results using an alternative nu- repr esentations of the free-space propagators, the reflec- 

merical procedure. Moreover, the method developed be- ^ mcthod and Schwinger > s spectral z _ forms g 

low gives also insights into the order of magnitude of the g|0|i|[n||J|]|||]. Section III contains a derivation 

coefficients of subleading divergences and finite parts. In of ^ expressions for thc 00-component of the canoni- 

principle, these can be calculated to any given accuracy, c&1 energy _ momeiltum ten sor, integrated over the angles, 



1 The free-space global divergence which renormalizes the for Klein-Gordon and Dirac fields satisfying the linear 
cosmological constant is, in this case, already subtracted. This boundary conditions of the M.I.T. bag model on a static, 
corresponds to the condition that the vacuum expectation val- spherical surface. Using the cavity mode representation 



ues vanish in the case of a free-space theory. of the propagators [115,23,24], the numerical evaluation is 
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done by summing over the cavity mode quantum numbers Smit,<i(x, y) — S°(x — y) — ij 3 S°(x — y T ) , (6) 



first and integrating the resulting spectral z-form subse- yT _ yi y2 _y3^ 

quently. An integration over a regularization volume yields 
the canonical regularized energy. Thereby an expansion in 
terms of the regularization parameter is derived from the 
expansion of the density into a Laurent type series . In 
section IV, the expression for the fermionic bag constant 
is derived and evaluated in analogy to section III. In the 
last section, we summarize our results and compare them 
with those available in the literature, and subsequently, we /qOO 
discuss problems and speculate on their possible solutions. 



where S° and A stand for the free-space propagators. 
The structure of the canonical energy-momentum tensor 
9^ u (sec for example ^5|) requires that, in order to obtain 
the (diverging) vacuum energy densities, some bilinear op- 
erators must be applied to the corresponding propagators. 
For a massless Klein-Gordon field, we obtain 



(x)) = lim d x od y o A{x,y) , (7) 
y—>x 



2 Vacuum energy in half-space field theories 



whereas, in the case of a massless Dirac field, the canonical 
vacuum energy density reads 



<£P(x)) = lim-Tr[(d x o-cV) 7 °] S(x,y) , (8) 

In this section a method for calculating the vacuum ex- y ~* x 1 

pectation value of the canonical energy density 6» 00 is illus- where ^(^ y ) and g( X) y ) denote the field propagators, 

trated for half-space field theories. Let us consider a mass- respectively 

less scalar field subject to Dirichlet (D) and Neumann (N) Let us now give a gener al illustration of the z-form 

boundary conditions, and a massless Dirac field satisfying method for the calculation of the vacuum expectation 

the linear boundary condition of the M.I.T. bag model yalue of local fidd bilinears . As an example, we concen- 

(MIT,q) Q trate on the expectation value of a Lorentz scalar in a 



scalar field theory, characterized by the bilinear This 
scalar quantity can be expressed in terms of the identity 



D: <i>{x)\ xe s = (1) 
N: d„n^{x)\ xeS = (2) 
MIT, q : (in^y - l)^(x)| Be s = , (3) {B(<f>(x),tf,{x))) = lim B x , y A a (x - y) . (9) 

n:= (0,0,0,1), S := {x\ x 3 = 0}. 

Here A (x — y) stands for the free-space Feynman prop- 
Using the reflection method the half-space propaga- agator, and B x>y denotes the point splitted version of 



tors read 2 The Dirac case can be treated in straightforward analogy. 

3 A contribution of a product of field derivatives 

A D (x, y) = A°(x - y) - A°(x - y T ) (4) d^ix) d„<p{x) to B, e.g., implies a contribution of c^cV to 

A N {x,y) = A a (x-y) + A°(x-y T ) (5) *W 
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In momentum space, the propagator is given by tions, we obtain 



1 /fl00/_ 3\\ _ 1 f j±„ (P ) n -2ix 3 p 3 

4°(p)oci, (10) (° ( x >)d--J2^J d P^TM G 



1 r d , p (p°r 



dz e~^ 2 /* z 



3 



and the application of B results in 

(2tt) 4 J e p 2 
= — f 

e -ip(x-y) h t 2\ . 32 7T 2 7n 

2 oc e"^*-") , (11) i J ° 

where b is a polynomial in n 2 (in a scale invariant theory b , .. . . . ... . 

while tor Neumann boundary conditions we have 

is just a power in p 2 corresponding to the mass dimension 

(e 00 (x s )) N = - 1 . (16) 

of B((f)(x), <j>{x)y). We now take the limit y — > a; and rotate n [x ) 

to Euclidean momentum space In the case of a massless Dirac field, subject to the linear 

boundary condition of the M.I.T. bag model for the quark 
P ° P ° field (see Eqs. we obtain 



The (Euclidean) denominator can then be elevated into ^6' 00 (x 3 )) =0 (17) 

This result is a consequence of the tracelessness of the 



' MIT,q 

an exponent 

1 

P z Jo energy-momentum tensor for this field 

and the remaining task is to integrate over the (Euclidean) 



dze~ zp , (12) 



momentum which involves (modified) Gaussian integrals 3 Vacuum energy in cavity field theories 

/dp 
J2n) • • •} e ^ = In this section, we shall calculate the canonical energy 

{-^ J_ I 1 „ _ q g ^3) density of the vacua of massless scalar and Dirac 

Calculating (9 00 ) for a massless scalar field, distorted by 



a Dirichlet plate at x 3 = 0, is now straightforward. The 
full propagator for this problem reads 



fields, 

confined to a static and spherical cavity with radius R. 
After a brief description of the methods, the regularized 
vacuum energy is calculated for each field. 



A(x,y) = A°(x -y)- A°(x - y T ), (14) 31 Massless scalar fields 

y T = (y° , y 1 , y 2 , — y 3 ) . In order to perform the free-space subtraction in the spher- 

ical symmetric case, we need to express the free scalar 

Renormalizing the cosmological constant, by applying the 

propagator in terms of angular momentum eigenstates us- 

bilinear operator of Eq. ([?]) to the difference of the half- 

ing the Rayleigh expansion for plane waves 

and free-space propagators, yields finite results except at 

,■ r p-ip(x-y) 

A°(r — v) = — — — I d 4 n - 

the boundary. In the case of Dirichlet boundary condi- v ' y > (27r) 4 J p 2 + iO 
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The result of the free-space subtracted canonical vac- 
uum energy density in the cavity is then 



— (4tt) z - — -j / (Fp , 
1 ; (2tt) 4 J y p 2 +i0 

x E i\-if 3i{w)ji>{pr') Y l>m (f) (18) 



Applying the bilinear operator of Eq. (0) to this represen- _ 1 f°° , 1 

~ 4 7rV2 J ^ Z 3/2 2^ 



tation of the propagator, and performing in turn the an- 
gular integration, the summation over m and a Euclidean 



rotation yields, after introducing the z-integration accord- _ 2(2/ + 1) 



ing to Eq. ©, 



dk k 2 (Jl(kr)) 2 e" zfc2 



4 TrVa J 



-i /-OO /"OO U / 



dfc k 2 {j t {kr)) 2 e~ zk " } . (22) 



x W) 2 ip- ( 19 ) K 2 r > 

io 

For the cavity part, we obtain the vacuum expectation 

value of the canonical energy density by applying the bi- Here the a^f jx) denote the scalar cavity modes for either 



linear operator of Eq. @) to the cavity mode representa- Dirichlet or Neumann boundary conditions, and the Af; 



D,N 



tion of the propagator (see Appendix |A.l|). A z-integration stan d for their normalization constants fB|,[||, respec- 

is introduced, which originates from a shift of the (Eu- tively, as explained in Appendix |A~l|. Using the plane- wave 

clidean) momentum squared denominator of the propa- representation of the free-space propagator implies a z" 3 

gator into an exponential, as discussed in the preceding divergence in the free-space part of Eq. Since finite 

section. Here the meaning of the term momentum differs sums over linearly independent functions cannot change 

somewhat from that of free space due to the boundary tne divergence structure common to all of these terms, a 

conditions, i.e. the analogue to the expression numerical evaluation of (# 00 (r)>D,N would diverge. 



P 2 = (p°) 2 ~ (p) 2 (20) The calculation reveals a substantial difference between 

the Dirichlet and the Neumann case. The Dirichlet bound- 
ary condition yields a z-form that is integrable at z = 0, 
whereas the Neumann boundary condition leads to a non- 
where w denotes the arbitrary (off-shell) energy, and s n ,i integrable z -3 / 2 divergence. We are able to show that this 



in free space, is in the cavity 

( Pn ,i) 2 = u; 2 - (e nA ) 2 , (21) 



stands for the energy of the mode, labelled by the ra- divergence is a global one, i.e. independent on r' := r/R. 
dial quantum number n and angular momentum quantum It therefore resembles another volume divergence (as does 
number I. the free-space global divergence), and hence it can be 



3.2 Massless Dirac fields 

Using Eq. (||), the mode summation representing the cav- 
ity Dirac propagator Jl^,^J|, the spherical representation 
of the free-space propagator, and introducing the Schwin- 
ger z-integral, we obtain, after an angular integration, 



300 



(r)) :=4tt ((9 00 (r)) 

MIT.q X ' 



dz ■ 



MIT,q 



2 TrVa J Q z 3 / 2 
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omitted through a renormalization of the cosmological 3.3 Numerical evaluation of (9 00 ) 

constant. 

Equations and (|2^) are suitable for a numerical eval- 
uation. The k integration is performed in the range from 
zero to £ max , the maximal energy eigenvalue used in the 
sum over cavity modes. In our computations it is typically 
of the size 200 /R. 

The calculation of (0 00 ) is done in two steps. At first, 
we compute the z-form for the corresponding field and 
linear boundary condition at a number M of points r' 
(M w 500 and r' := r/R). It thereby proves convenient to 
make the variable substitution z — y 2 resulting in a pure 
y~ 2 divergence in the y-form. In a second step, we inte- 
grate the regular part of the y-form and hence determine 
> the energy density as a function of the position r' . 

In Figs. [I] and |2|, the y-forms for (8 00 ) are displayed 
for the scalar Dirichlet and the Dirac case, respectively. 
It is shown numerically that the region of small y, where 
the form is practically zero, moves to the left with in- 

(23) creasing maximal energy £ max , whereas points to the right 
of this region remain unchanged. The y-form for (9 00 ) in 
the scalar Neumann case contains a (under an increase 

(24) of e max ) stable and global, i.e. independent of r' , y~ 2 - 

(25) divergence (see Fig. ||). The subtraction of this divergence 
amounts to renormalizing the cosmological constant and 
results in a y-form shown in Fig. [| 

The error /(r',e max ) of (9 m ) due to the truncation of 



E^E / ^ 1^)n,nA e" Ze " 

, k n, fi 



, 7T J 



dk k 2 {ji{kr)f e 



dz ■ 



2 TrVa J Q z 3 / 2 

+ (ir(k«,K|r)) 2 )e- ; 

-E - ( 2Z + 1 ) I dkk 2 { 3l {kr)) 2 e~* k2 



where 



J := | K | 
I ■= \J\ 



sgn k 



1 = 1 — sgn k . 



(26) 



Here, q n ,K,n{ x ) denotes the Dirac cavity mode, labelled 

with the radial quantum number n, the Dirac quantum - A • ^ /looh j /Roh u 

H ' H the sum over cavity energies m Eqs. (22) and (E3J) can be 



number k, and the angular momentum projection fi, as determined by varying thc cut _ off cnergy £ ^ We can 

safely estimate the upper bound for /(r',e max ) at 10~ 6 , 



discussed in Appendix A. 2. 



valid for all r' and all fields and boundary conditions under 
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consideration. As an example, Fig. |i]shows {6 00 (r')} for a The numerical errors of the coefficients c_4, . . . , c_i can 
Dirichlct scalar field. 

Following Deutsch and Candelas |1, we may expand 
(9 m ) into a Laurent type series around r' = 1 



be estimated by varying the interval 
<5Ln <S'< 5' 



(6 00 }(5') 



1 



c_ 4 (5')- 4 + c-3 (ST 3 + --- 
+ C _! (ST 1 +c + -.. + c N (r') N 

N > , 8' := 1 - r' . 



,(27) 



To extract the coefficients of the negative powers in Eq. 
@), the function 



In our computations we vary S' min from 0.05 to 0.15 and 
8' max from 0.1 to 0.2 to obtain errors of about 1%, 10% 
and up to 130% for the coefficients of the leading, next to 
leading, and the weakest divergences, respectively. 

In the literature, there have been several attempts to 
extract boundary divergences of (9 oa ) or the canonical 
energy E(e') using the asymptotic properties of analytic 
functions |^[l(].|l2|,^6]]. The drawback of these expansions 
lies in the fact that one does not know the errors of the 
analysis In our calculation of the regularized energy, 
we start with a determination of (6 m ), where the only 

source of substantial errors is the fitting procedure to the 

tion of the positive power coefficients in Eq. (E7|), wc fit T , , TT , 

1 1 H v l — V Laurent type series. However, these errors can be esti- 



£'\-4 



E^{8') := c_4 (5' 
is fitted to the calculated curve in the interval 



min — — max ' 



where S' ; „ and 5' av are close to zero. For the determina 

Illlll IllciX 



the polynomial 



E (/) :=«(, + «! r' + ... + (/) 



mated and made smaller in more extensive numerical cal- 
culations. 



to the calculated curve within the interval 
< r' < r' -=1-8' 



3.4 The regularized canonical vacuum energy 



Following Bender and Hays pTj| , we regularize the canon- 

A value of N = 9 is used to obtain negligible fitting errors. . n ^ , . , ,. //inn/ iw i 

00 ° ical energy E by integrating \0 m (r )) only to an upper 



Fitting E L (S'(r')) to the polynomial 



limit of 



b + bi r + ■ ■ ■ + b N (/) 



(l-s'), e':-- 



R 



. . , r inoo\ j where e denotes the distance from the boundary. Using 

within the above interval, the expression for (6' uu ) reads J ° 

Eq. (p7|), we obtain 



;0 oo (r')> - E L (S'(r')) + AE(r') 1 
AE(r') := co + Cl r' + • • • + c N (r') N , (28) 
d := at -h, (0 < i < N) . 



E(e') 



- j £ dr' r' 2 (9°°(r')) 

R Jo 
l 

R 



I -3 . ~ / -2 _ , -1 

c_3 e + c-2 £ + c_i e 
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ergy ||]. In bag model calculations B is determined using 
the minimization condition 
d 



+ ci og log e +c n + 0{e') 



(29) 



The coefficients {c} of Eq. (p7|) and the coefficients {c} of 



Eq. (|2g) are related by dR E ^ R > ° 



5-3 = | C_ 4 , 



where E(R) denotes the total energy of the bag. Milton 



5-2 = -C-4 + 2 c -3, 



i c [[12| suggests that one should calculate the bag constant 

(30) 

g _ c 4 _2c 3 + c 2 from first principles, by interpreting it as a vacuum ex- 



gj _ _ c ^2 c 2 _ c 1 pectation value. The nonlinear boundary condition of the 

r , rc ■ , r ,i model then serves as a definition for B. Here we will carry 

for the divergent terms, whereas the coefficients of the J 

r / j ii , i i out the calculation of this vacuum expectation value for a 

positive powers of e depend also on tne truncation number ^ 

AT c . i • , , , • ,i i- ■ / n massless fernfion field, subject to the boundary conditions 

N. Since we are only interested m the limit e 0, the ' J J 

i uj-i-i a;-i-cr~i-\m-~ • u of the M.I.T. bag model on a static sphere with radius R. 
only substantial coefficient of {Cj \ j > 0) is Co, given by ° r 

jv The expression for i?MiT,q reads 



1 3 3 1 

For the fermionic case, we may assume c_4 = 0. This 
is so for fermionic quantum fields because the canoni- 



B M iT, q := ~{d r (ipip)) 

1 r=R 

and implies that the differential operator 



(32) 



cal energy- momentum tensor coincides with the improved — lim — Tr (<9r x +<9 ry ) (33) 



tensor, and hence it is traceless [||. 

Table |^ contains a list of the coefficients {c} for the 
various fields and boundary conditions, where the errors 
are determined from the errors of the coefficients {c} ac- 
cording to the rules of error propagation. Note that the 
leading divergences for the scalar Dirichlet and Neumann 
fields exhibit a behavior analogous to that of the half-space 
problem, i.e., they are equal in magnitude and carry op- 
posite signs. 



y^x 2 

should be applied to the difference of the cavity and the 
free-space propagator. We expect i?MiT, q to be infinite. 
In order to regularize it, we compute it at some interior 
point a distance e away from the boundary, rather than 
at some point on the boundary, as Eq. (|32|) demands. For 
the angular integrated version i?MiT,q(£') of i?MiT,q(£')i 
we arrive at 



S M iT,q(e') := 47r B M iT, q (e') 

1 f°° 1 

= -zip / dz 7m E( 2J + !) E 

4 The fermionic bag constant x | J ' ( ^j r) (j ji^(\e n , K \r) - (I + 1) ji+i(\e n , K \r)) 

jj(\s n K \r) - - 1 

The bag constant B is introduced into the Lagrangian to 21 + 1 ('^- 1 (l en ' K l r ) — ^ ^)^+i(l £ n,Kl r )) r ; 

achieve conservation of a Poincare generator, i.e. the en- (34) 
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with inition via the nonlinear boundary condition was carried 

s = — , r := 1 — e . out. Numerical results for the densities have been fitted 
R 

rpi • r i • /fTlx • rp u — n c to Laurent type series in powers of the distance to the 

I here is no free-space part in hq. ( |34| ) since ir 7^ = for J v 1 



H = 0, . ..,3. 



boundary. The expressions obtained in this manner could 



The result of the calculation of the coefficients {c} be inte S rated within a regulation volume to yield an 

, , r ,i_ , , 1 , expansion of the energies in terms of the regularization 

m the expansion ot the regularized bag constant energy ± e> 

n R / /\ .,, parameter. Numerical errors of the coefficients in these 

E MIT, q ( £ ) Wlth 

expansions turn out to be less than 1% for the leading 



R 3 

^MIT,q( e ') : = — (1 - e ') x B M lT,q(s') 



"* singularity and up to 100% for the weakest divergences. 

1 



R 



, -4 - , -3 . _ , -2 

C_4 £ + C_ 3 e + C_2 E 



In general, the finite part Co could be determined to about 

c_i e' 1 + So + O(e') (35) 50% accuracy, using this method and the stated compu- 

_ tational effort, 

is displayed in Table 0. Here we also list the set of co- 

We can compare our results directly with those in the 

efficients in the expansion of the total vacuum energy 

, . , literature. Bender and Hays [fl(J| have calculated the lead- 

E&rje') given by ^ 

l 1 ing divergence of the canonical part of the vacuum energy 

£M°iV q (£') == / dr'(rr{(0 OO (r'))MiT, q -B M1T . q (e')} . . . 

,4 Jo I ' ™ J using a Green s function method. They introduced their 

— ^MIT,q(e ) — E MIT q (e ) . (36) regularization in the same fashion as we do, but they re- 

lied on analytical formulae for the asymptotic expansion 
5 Summary and Discussion of Bessel functions. For example in the fermionic and the 

Dirichlct scalar case, the comparison is as follows: 

The main purpose of this paper was to investigate the 

effect of a static and spherical boundary on the canon- Bender & Hays : c_ 2 ,q = —j^n ~ — 0.0027 , 
ical vacuum energy densities of otherwise free massless 



= ~ -0.0133 



Klein-Gordon and Dirac fields. Thereby the linear bound- this work : c M 2 T ' q = -0.01058 ± 0.45% 



(37) 



ary conditions of the M.I.T. bag model H have been used. ^-3 = 0.01324 ± 0.1% 

A Green's function method, that is based on the eigen- There are, of course, disagreements, but the sign error of 

mode representation of the propagator and the Schwinger Bender and Hays' result for the scalar case has already 

parametrization of the Euclidean momentum squared de- indirectly been pointed out by Milton [p^[ . In his work 

nominator, has been used to obtain numerical results for he mentions an overall sign error in the vector field mode 

the vacuum energy densities. For the fermionic field, a cal- sum of Ref. @. Since the mode sum of the transverse 

culation of the regularized bag constant based on its def- electric vector field is up to the I = contribution equal 
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to that of the scalar Dirichlet case (compare Eqs. (2.27) The small error of the leading infinite terms allows 
and (2.15) in Ref. fuj), an overall sign error implies a sign us to compare the results for the canonical and the bag 
error in the result of the scalar Dirichlet case. The factor constant part in the fermionic case. A cancellation of the 
of four difference in the fermionic case might be due to an infinities between these two contributions to the total en- 
omission of the trace of 7q as required by Eq. (^). ergy does not occur, since there is a nonvanishing quartic 

m 1^11 rFT^i i r~i i r divergence in the bag constant contribution, whereas the 

Olaussen and Ravndal |E6| used a Green s function ° ° 

,11, i n , , , ,. i leading infinite term in the canonical part is quadratic 

method to calculate the electromagnetic canonical vacuum ° r -i 

, ., r i • ,,, , , U1 (see Table 0). Moreover, the coefficient en in EEt T _ is 

energy density for a spherical bag, and we are thus not able \ a/ > u Mii,q 

to compare their results with ours. quite large. 

There are two common types of regularization proce- One may speculate that, in a gauge theory calculated 
dures. For example Milton []T2|,fl3| calculated the canonical perturbatively, infinities arising from the canonical and 
fermionic and vector field vacuum energy using a tempo- the bag constant part of the energy should cancel each 
ral regularization, whereas Bender and Hays |lC| ] worked other, leaving a meaningful finite part. We plan to in- 
with the same volume regularization as we did. Tempo- vestigate this matter in the near future. If there is still 
ral regularization derives from the idea of point splitting no hint for a cancellation, one would argue that the in- 
as a means of cutting off ultraviolett contributions to the troduction of a constant bag energy density B is a too 
energy density at each point in the cavity. On the other naive device to achieve the conservation of the Poincare 
hand, the volume regularization used above avoids the in- generator, i.e. the energy. Perhaps, a locally conserved 
tegration of boundary divergences of the energy density, energy-momentum tensor (if at all definable on physi- 
Therefore a strict comparison of results obtained in those cal grounds) would reveal finite vacuum energies. On the 
two different regularization schemes should not be made, other hand, using the canonical energy-momentum tensor 
One would, however, expect the same divergence structure together with soft boundaries, i.e. practically confining po- 
with different coefficients. An order of magnitude check of tentials (as for example a harmonic oscillator potential), 
Milton's result in comparison with our result for c_2 in could possibly bypass the global vacuum infinities. The 
the fermionic case |l2| reveals less than a factor of 10 dif- drawbacks of this method are the enormous numerical ef- 
ference. The large deviation in the finite parts (a factor of fort and the necessity of a cutoff parameter to ensure con- 
about 10 2 ) can very well be a consequence of the use of finement. Moreover, it is quite unsatisfactory, that this 
asymptotic expansions in Ref. |l2]|. In fact, as Milton him- additional parameter would have to be determined ex- 
self points out in the conclusions of Ref. Jl3| , even existing perimentally. Along the same lines, there have been sug- 
logarithmic divergences are then not recognized. gestions to introduce phenomenological parameters in the 
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expression for the bag energy. These parameters could ab- The energy eigenvalues e n ,i depend, of course, on the cho- 

sorb the divergences. However, they should be determined sen boundary condition. In the Dirichlet case, we have the 

experimentally which is again unsatisfactory. eigenvalue equation 

A possible extension of the work done here would be 
a more extensive numerical calculation (determination of 
the energy density at points closer to the boundary) to 
achieve higher precision for the coefficients of the sublead- 

ing divergences and the finite part. ljl(^n,lR) ~ e n jRji +1 {e n jR) = . (41) 



ji(e n ,iR) = , (40) 

and in the case of Neumann boundary conditions the eigen- 
value equation is 



The normalization constants Af^i ano - -^Cz are given as 



/w-d.n _ 
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(42) 



the Foundation for Fundamental Research (FFR) and the Gra- 7V D ^ N = (431 
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duiertenkolleg "Struktur und Wechselwirkung von Hadronen 

und Kernen" are gratefully acknowledged. A. 2 Massless fermion fields 



The cavity modes for massless fermion fields satisfying 
A Cavity modes the linear boundary condition of the M.I.T. bag model 

are given by the Dirac spinors 

Here, the eigenmodes and propagators for massless scalar 

q nK r x)= ( 9n, K (r)xW) \ {u) 
and Dirac fields in a static spherical cavity of radius R are \*/n,K( r )X-K;(^) / ' 

given |p3|, p4] , [2^1 . where Xk(^) i s ^ ne usual two-component spherical spinor. 

Here n, k, and /i denote the radial, Dirac, and magnetic 

A 1 Massless scalar fields quantum numbers respectively, and the radial functions 

9nA r ) and / n . K (r) are given in terms of the spherical 

The cavity modes of the massless scalar fields arc given as „ , r , . . , 
J ° IJessel functions ji by 

W D ' N Kf 
<Pn,lA^) = -j^jy i 3i(£n,ir)Y lfl (r) . (38) 9nA r ) = ^| hipn,^) (45) 

/„ f r ) = M,,ttSgn(rc) S gn( K ) _ 
A mode representation of the propagator reads as R?l 2 



A(x,x') = i ( x )#M )A t( x ') 



The discrete momenta p n ^ K in Eqs. (45) and (46) are de 



n : i tll termined by the linear boundary condition 

•<&j e -M*o-*'o) 

2^oj 2 - {e n ,i) 2 + i0 ' (39) (h-r + l)q n ,KA r )\r=R = (47) 
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which leads to the eigenvalue equation 8. M. Bordag, E. Elizalde, and K. Kirsten, J. Math. Phys. 

37, 895 (1996). 

jl{xn, K ) +sgn(n)sgn(/t) j r (x„, K ) = . (48) 

9. D. Deutsch and P. Candelas, Phys. Rev. D 20, 3063 (1979). 
The normalization constant Af n>K in Eqs. @ and (|(| is in. C. M. Bender and P. Hays, Phys. Rev. D 14, 2622 (1976). 
given by 11. k. A. Milton, L. DeRaad, and J. Schwinger, Ann. Phys. 

115, 388 (1978). 



ji (x n ) 



(49) 



12. K. A. Milton, Phys. Rev. D 22, 1444 (1980). 
Here we have introduced the dimensionless energy and 13 R A mto ^ phyg Rev D 22 M41 (lgg0) 

momentum parameters u T H Hangson and R L Jaffij) phys Rev D 28 g82 



(1983). 

15. A. J. Stoddart, Ph.D. thesis, University of Cape Town, 
1989. 

The cavity propagator for massless fermions can be rep- 16. G. U. Schreiber and R. D. Viollier, Ann. Phys. 215, 277 



X n ,n — Pn.nR-i (50) 

uj UtK = sgn (n) x„ iK . (51) 



resented in terms of these cavity modes as (1992). 

duj e -iu,(x -x' ) 17. G. U. Schreiber and R. D. Viollier, Phys. Lett. B 279, 131 



KVfj, J ™ (1992). 



S (x,x') = iy^gn, K ,At( x )gn,«,M( x/ ) J 

(52) 

where the usual Feynman prescription for the poles is 
used. 
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Fig. 1. y- forms for (8 00 ) at three different locations r' = r/R 
for a massless scalar field fulfilling Dirichlet boundary condi- 
tions on a static sphere with radius R. 
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Fig. 2. Massless fermion field y-forms at three different lo- 
cations r' — r/R. The field modes fulfill the linear boundary 
condition of the M.l.T. bag model on a static sphere with ra- 
dius R. 
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Table 1. The coefficients {£} for the divergent and the finite parts of the canonical vacuum energy E(e') for massless scalar 
Dirichlet (D) and Neumann (N), and Dirac (MIT,q) fields. 



C-3 



C-2 



C-l 



Clog 



Co 



D 0.0132435 ± 0.1% 
N -0.01318 ± 0.5% 

MIT,q 



-0.0232 ± 10% 
-0.0012 ± 60% 
-0.01058 ± 0.45% 



-0.0065 ± 40% 
-0.085 ± 40% 
-0.35 ± 100% 



-0.0052 ± 10% 
0.64 ± 15% 
-0.73 ± 100% 



-10.2 ± 50% 
45 ± 60% 
0.16 ± 20% 



Table 2. The coefficients {£} for the divergent and the finite part of i?MiT,q an d the total vacuum energy i?Mrr,q(£')- 





C-4 


C-3 


£-2 


C-l 


Clog 


CO 


-C'MIT.q 


-0.1593±0.2% 


0.375±1% 


-0.119±5% 


0.11±20% 





-7.2xl0 5 ±50% 


j-itot 
^MIT.q 


0.1593±0.2% 


-0.375±1% 


0.108±6% 


-0.46±85% 


-0.73±100% 


7.2xl0 5 ±50% 
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Fig. 3. Unsubtracted and subtracted y- forms for (6 00 ) at Fig. 4. Vacuum energy density for a massless scalar field ful- 

r' — 0.5 for a massless scalar field fulfilling Neumann bound- filling Dirichlet boundary conditions on a static sphere with 

ary conditions on a static sphere with radius R. The subtracted radius R. 
divergence is of the form m/y 2 , where m = 0.8463. 



